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Abstract 

In recent papers, we have studied the quantum-mechanical decay of a Schwarzschild-like 
black hole, formed by gravitational collapse, into almost-flat space-time and weak radiation 
at a very late time. In this recent work, we have been concerned with evaluating quantum 
amplitudes (not just probabilities) for transitions from initial to final states. In a general 
asymptotically- flat context, one may specify a quantum amplitude by posing boundary 
data on (say) an initial space-like hypersurface Ej and a final space-like hypersurface E^ • 
To complete the specification, one must also give the Lorentzian proper-time interval be- 
tween the two boundary surfaces, as measured near spatial infinity. We have assumed that 
the Lagrangian contains Einstein gravity coupled to a massless scalar field <fi , plus pos- 
sible additional fields; there is taken to be a 'background' spherically-symmetric solution 
(7^,$) of the classical Einstein/scalar field equations. For bosonic fields, the gravita- 
tional and scalar boundary data can be taken to be gij and <fi on the two hyper surfaces, 
where gij = 1,2,3) gives the intrinsic 3-metric on the boundary, and the 4-metric is 

9[iv (a*) v = 0,1,2,3), the boundary being taken locally in the form {x° = const.}. The 
classical boundary-value problem, corresponding to the calculation of this quantum am- 
plitude, is badly posed, being a boundary-value problem for a wave-like (hyperbolic) set 
of equations. Following Feynman's +ie prescription, one makes the problem well-posed 
by rotating the asymptotic time-interval T into the complex: T — > \T\ exp(— i9), with 

< 9 < 7r/2 . After calculating the amplitude for 9 > , one then takes the 'Lorentzian 
limit' 9 — > 0+ . Such quantum amplitudes have been calculated for weak s = (scalar), 

s = 1 (photon) and s = 2 (graviton) anisotropic final data, propagating on the ap- 
proximately Vaidya-like background geometry, in the region containing radially-outgoing 
black-hole radiation. In this paper, we treat quantum amplitudes for the case of fermionic 
massless spin-| (neutrino) final boundary data. Making use of boundary conditions orig- 
inally developed for local supersymmetry, we find that this fermionic case can be treated 
in a way which parallels the bosonic case. In particular, we calculate the classical action 
as a functional of the fermionic data on the late-time surface E p ; the quantum amplitude 
follows straightforwardly from this. 

1. Introduction 

This paper treats quantum amplitudes involving neutrinos (s = |) in the context of 
nearly-spherical Einstein/massless-scalar gravitational collapse to a black hole, following 
the approach of [1-9]. We first summarise the main ideas at the bosonic level. Writing 

9fj,v (fJ-iV = 0,1,2,3) for the 4-metric and hij = gtj = 1,2,3) for the intrinsic Rie- 
mannian spatial 3-metric on a surface {x° = const.}, one poses appropriate boundary data 

(h%j,<f>)i,F ( sa y) on an initial space-like hypersurface Ej at x° = and on a final surface 



1 



Yi F . For simplicity, we assume that both Ej and Hp are diffeomorphic to IR 3 . The large 
(real) Lorentzian time-separation T between the surfaces is measured at spatial infinity. 
Again, for simplicity, the initial data are taken to be exactly spherically sym- 

metric, corresponding to a diffuse, slowly-moving, initially-spherical configuration. The 
space-time geometry near the final surface E^ is taken to be approximately a Vaidya 
metric [4,5], describing the region which contains radially-outgoing radiation due to the 
collapse to a black hole, except for small perturbations in the gravitational and scalar data. 

We assume that there is a 'background' Lorentzian-signature spherically-symmetric 
solution of the coupled field equations, with 4-metric 7 Mi , and scalar field $ . If the data 
for the gravitational and scalar fields are perturbed slightly from spherical symmetry, the 
resulting non-spherically-symmetric solution can be written in the form 

= 1^ + eh$ + e 2 h$ + ... , = $ + e^ 1 ) + e 2 <j>V) + ... , (1.1) 

where e is a measure of the size of the perturbation. We then ask (in this bosonic case) for 
the quantum amplitude to go from asymptotically-flat initial data (gij , 4>)j on the initial 
hypersurface E / to corresponding final data (g^ , (j))p on the final hypersurface E p . For 
a complete specification of the boundary data involved in the quantum amplitude, one 
also needs the (Lorentzian) time interval T between E/ and E_p , as measured at spatial 
infinity. 

Naively, one might expect there to be a simple (semi-)classical analogue of this quan- 
tum amplitude, as specified by the above boundary data. This would involve finding a 
classical solution of the coupled Einstein/scalar field equations, which agrees with the 
boundary data, and such that the Lorentzian proper-time separation at spatial infinity is 
indeed T . Unfortunately, this will only be possible for very special choices of boundary 
data, having measure zero among all possible boundary data. Equivalently, as is well 
known [2,10,11], the boundary- value problem for hyperbolic (wave-like) partial differential 
equations is not well posed. 

The procedure for understanding our (Lorentzian) quantum amplitude semi-classically 
is a little more complicated. Following Feynman's +ie proposal [12] , we rotate the Lorentzian| 
time-interval T into the complex: T — > |T| exp(— i6), with < 9 < n/2. In the extreme 
case 9 = n/2, the Riemannian region with positive-definite 4-metric is reached. In 
the Riemannian case, the field equations are 'elliptic modulo gauge', and one might reason- 
ably expect good existence, uniqueness and analyticity properties for classical solutions, 
under suitable conditions on the boundary data. In the intermediate complex case, with 
< 9 < 7r/2 , the interval T and any classical solution will involve the complex num- 
bers non-trivially. In such 'complexified elliptic' cases, it frequently turns out that the 
boundary- value problem is strongly elliptic, up to gauge, as defined and discussed in [13]. 
Strongly elliptic partial differential equations are again liable, as in the real elliptic case, 
to have good existence, uniqueness and analyticity properties [13]. There is reasonable 
ground for expecting that our classical boundary-value problem will have an (analytic) 
complex solution, for < 9 < n/2, at least under suitable conditions on the boundary 
data. 

One can then form the (complex) classical action S c \ a ss , for the given real boundary 
data [here (g^ , </>)], as a function of 9 . Similarly, provided that the quantum field theory 
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under consideration has meaningful quantum amplitudes (see below), there will be one- 
loop, two-loop, . . . factors Aq , A\ , Ai , . . . , again depending on the boundary data and 
on 9 . That is, the quantum amplitude will have an asymptotic semi-classical expansion 
of the form 

Amplitude ~ (A + HA 1 + ti 2 A 2 + ... ) exp(zS c i ass /fr) . (1.2) 

Provided that the frequencies involved in the prescribed boundary data are well below the 
Planck scale, one expects to be able to approximate the quantum amplitude by (const.) x 
exp(z5'ciass) , in units such that h = 1 . For dimensional reasons, any loop corrections will 
then be extremely small, compared with the leading term Aq exp(i5 , c i aS s) • Finally, one 
recovers the Lorentzian quantum amplitude, which is defined to be the limit as 6 — > + 
of the ^-dependent quantum amplitude. Provided that one is away from the Planck scale, 
then it is sufficient to compute limg^o + exp(i$ciass) • 

The only field theories which can give meaningful quantum amplitudes in our case, 
with the property (1.2), are locally supersymmetric, that is, supergravity theories with 
or without supermatter [11,14-16]. The simplest theory of N = 1 supergravity coupled 
to supermatter, which contains Einstein/massless-scalar theory within its bosonic part, is 
given in [17]. For this theory, the scalar field has become complex, with massless spin-i 
partner, and of course there is also a spin-| gravitino. It should be understood that, in 
the discussion above concerning quantum amplitudes for purely bosonic boundary data, 
the quantum amplitudes refer in fact to the full locally-supersymmetric theory. 

It is then natural to ask for quantum amplitudes involving boundary data which 
include non-trivial spin-i or spin-| contributions. In this paper, we shall compute the 
amplitude for non-zero spin-| data given on Ep, assuming that the spin--| data on are 
zero. In particle language, this would refer to the production of neutrino-antineutrino pairs. 
The resulting classical fermionic fields and action will then be elements of a Grassmann 
algebra, as usual in the holomorphic representation for fermions [11,18-20]. (In Lorentzian 
signature, this point of view has also been taken in [21], concerning the Cauchy problem 
for N = 1 supergravity.) In the present paper, the spin--| field is described in terms of 
2-component spinor fields [22,23], the unprimed spinor field 4> A (x) and the independent 
primed spinor field 4> A (x) , not related to 4> A (x) by any conjugation operator. The fields 
4> A (x) and (p A (x) are taken to be odd (anti-commuting) elements of a Grassmann algebra 
[18-20] - see Sec. 2 below. As discussed further in Sec. 3, this is the spin--| description 
appropriate to Euclidean signature (Riemannian geometry), and it is also the description 
appropriate to the 'Riemannian rotated-into-the-complex' case here, with < 9 < 7r/2 , 
and corresponding to a well-posed classical boundary-value formulation. The fermionic 
action (3.1) below is quadratic in spin-| fields or their derivatives, and so is an even element 
of the Grassmann algebra. This gives a 'field description' for the massless spin-| theory in 
our case, which allows a treatment of the fermionic quantum amplitude which is roughly 
in parallel with the previous treatments for bosonic fields [5,6,8,9]. The holomorphic 
representation being used here for fermions can straightforwardly be related to a (possibly 
more familiar) description in terms of particle states, as is outlined in Sec. 2. 

The basic fields and the spin-| action are discussed in Sec.3. The action has the form 
Si = Si v + Si B , where Si v is the standard volume expression, and Si B is a boundary 
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contribution, which depends on the particular boundary conditions being adopted. Natural 
local boundary conditions are taken, related to local supersymmetry [24-26]. These again 
allow a parallel between our s = | treatment and that for other (bosonic) spins. At 
a solution of the classical massless s = \ field equations (Weyl equations), the volume 
contribution Si v vanishes, whence Si c i &ss = Si Bc i ass . 

Some basic properties of the Vaidya space-time are summarised in Sec. 4. After the 
initial gravitational collapse, a large part of the classical Einstein/massless-scalar space- 
time can be approximated by a metric of Vaidya type [4,7,27], with its mass decreasing 
extremely slowly with time. This describes the region of the space-time containing the 
emitted radiation. Almost all the spin-i (or other spin) modes of physical interest will 
evolve adiabatically in this background. 

In Sec. 5, we review the separation of variables for massless spin-i fields in a Schwarzschild-| 
like background. This leads to a treatment of the classical action Si Bclasa in which local 
boundary data are specified on the final surface T,p , given, for simplicity (as above), that 
trivial fermionic data are taken on Ej . The form of the resulting classical action, as a 
functional of the final boundary data, is compared with that for massless spin-0 perturba- 
tions, found in [5,6]. The form of the spin-| Lorentzian quantum amplitude can then be 
deduced straightforwardly from the definition lim#^o+ exp(i5 , c i aS s) • A brief conclusion is 
contained in Sec. 6. 

2. Holomorphic representation for fermions 

Suppose, at first, that one had a Fermi system with only one degree of freedom [18,19], 
rather than a quantum field theory. In the holomorphic representation [18-20], one would 
begin by considering two variables a* and a , which are odd (anti-commuting) elements of 
a Grassmann algebra: 

a*a + aa*=0; (a*) 2 = ; a 2 = . (2.1) 

A typical element of the resulting Grassmann algebra (over the complex numbers) has the 
form 

/(a*, a) = f 00 + /oi a + / 10 a* + f n aa* , (2.2) 
where /oo > foi ■> fio an d fn are complex numbers. Functions of the form 

/(a*) = fo + fia* , (2.3) 

with f and fi being complex numbers, are called holomorphic functions; these describe 
the state vectors of the system and form a two-dimensional space. 

The operators a* and a, acting on state vectors given in the form (2.3), are defined 

by 

a* (/(a*)) = a* /(a*); a(/(a*)) = ^/(O , (2-4) 
where fermionic differentiation is defined by 

^ (/.+/-•) = h • (2.5) 
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Then, one can verify the standard anti-commutation relations 



a*a + aa* = 1; (a*) 2 = 0; (a) 2 = . (2.6) 
Berezin integration [18-20,28] is defined through the rules 

J a* da* = 1; J ada = 1 ; J da* = ; J da = . (2.7) 
The inner product of two state vectors is defined to be 

(/i,/ 2 ) = j{h{a*)Y f 2 {a*)e- a * a da*da , (2.8) 
This inner product is positive-definite. Indeed, the monomials 

V>o = 1; V>i = «* (2-9) 

are an orthogonal pair, each of unit length. One may interpret ipo as giving a no-particle 
state and ipi as giving a one-particle state; similarly, the operators a* and a give creation 
and annihilation operators. In this way, one can replicate the usual structures familiar in a 
particle description of fermions. But, in the holomorphic representation for fermions, the 
wave functions are encoded through the 'complex-analytic' form f(a*) = f + fi a* , which 
makes the parallel with bosonic systems much closer; in the holomorphic representation 
for bosons [18-20], wave-functions take the form f(a*) , where, in the bosonic context, a* 
is simply a (commuting) complex variable. 

This representation can be generalised to fermionic field theory [18,19], by considering 
a collection of Fermi oscillators interacting with an external field. In our case (see below), 
we fix half of the spinorial quantities (4> A (x) , <p A (x)) at each point x on £/ , and similarly 
half at each point on , in order to have a well-posed classical boundary- value problem 
for the massless Dirac (Weyl) equation. This amounts to choosing an initial (and a final) 
wave-functional, which is a product of delta-functions in the fermionic variables. In the 
simple Fermi example of Eqs.(2.1-9), one has (for example) d(a*) = a* . That is, in our 
case, the wave-functionals specifying the initial and final data are (formally) products, 
over all points x , of suitable functions of the type (2.3). This holomorphic approach 
for fermions, of course, exploits the connection between the quantum and the classical. 
Similarly for bosonic fields, where the Bargmann-Fock holomorphic representation is known 
as the theory of coherent states in quantum optics [18-20,28]. 

3. Massless spin-^ fields 

We are concerned with the quantum amplitude to go from initial data specified on an 
initial space-like hypersurface £/ , at time t = (say), measured at spatial infinity, to final 
data on a hypersurface Ej? , at a very late Lorentzian proper time T , again measured at 
spatial infinity. As described in Sec.l, we rotate the time-interval T into the complex: T — > 
\T\ exp(— z6>) , where < 9 < tv/2. In this case, the classical boundary- value problem is 
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expected to be strongly elliptic, up to gauge, leading to a complex classical solution (g^ , 0) 
for Einstein/scalar theory. From this, one can recover the Lorentzian quantum amplitude 
via Feynman's +ie procedure, by considering the semi-classical amplitude in the limit 
9 — > 0+ . As mentioned above, the Lorentzian-signature boundary- value problem (9 = 0) is 
badly posed [2,10,11], even for linear systems of equations - in general, there is no solution, 
or, if there is a solution, it may not be unique. It is for this reason that calculations in 
quantum field theory (along the lines just given) must proceed in the 'Euclidean direction', 
that is, by taking Im(T) < . 

For quantum purposes, then, we should think of the massless spin-| field from the 
Euclidean viewpoint. In the simplest case, one has an unprimed 2-component spinor field 
<p A and an independent primed spinor field <p A [11] (this reference describes the spinor 
conventions being used here). In addition to the familiar Einstein/massless-scalar action 
of Eq.(3.4) of [5], possibly including also a Maxwell contribution as in Eq.(2.1) of [8], the 
(Lorentzian) action 5 now contains an extra piece 

Si = - % - f d 4 x e j> A ' e AA r {D^ A ) + 'h.c' + Si B . (3.1) 

Here, as in Sec. 6 of [8], e AA ? denotes the spinor form of the tetrad e a M (a = 0, 1, 2, 3) , and 
we define e = det(e a M ) = (— g)z , where g = det(g fJ/ ^) . The spinor- valued connection 

one-forms u AB = u/" 45 ^ and their tilded counterparts u A ' B ' = uj^ A B \ are defined 
in Sec. 2. 9. 2 of [11]. The covariant derivative on spinors is defined, with the help of the 
distributive law, via 

D^ A = d^ A + u A Bll <P B , (3.2) 
D^ A ' = d^ A ' + u A B ,^ 4> B ' , (3.3) 

Also, the expression 'h.c' in Eq.(3.1) denotes an operation in which (f) A is replaced by 
<p A , and (j) A by 4> A , and in which the order of fermionic integrands is reversed. This 
is necessary, since fermionic quantities such as <p A or 4> A must be taken to be odd 
elements of a Grassmann algebra, anti-commuting among each other (obeying, for example, 

as in Sec. 2. The final element of Eq.(3.1) is a boundary contribution 
, which depends on the fermionic boundary conditions chosen (see below). 

On varying the independent quantities <p A , 4> A in the action (3.1), one obtains the 
Weyl (neutrino) field equations 

e AA r D^ A = , (3.4) 
e AA r D^ A ' = . (3.5) 

Next, consider possible boundary conditions. It will again be simplest and most symmet- 
rical to take local boundary conditions, as in the cases of s = [5,6], s = 1 [8] and s = 2 
[9]. By analogy with Sec. 6 of [8] and Sec. 10 of [9], define, at the boundaries E/ and : 

$ A (x) = eV2n A A/ ^ A '(x)+(l ) A (x) , (3.6) 
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where e = ± 1 and n AA is the spinor version of the unit time-like future-directed vector 
n M , normal to the surface Ej or Hp . If one regards the right-hand side of Eq.(3.6) as 
being given by a 4 x 4 matrix operator P e , acting on a column vector with components 
((f) , cf) 1 , 4>° , (j) 1 ) , one finds that P+ and P- are projection operators, obeying 

P+ + P- = 1 , = P_ P+ = , (P e ) 2 = P e . (3.7) 

The use of such local boundary conditions is suggested by the work of [24-26]. Given a 
choice of e, take the boundary contribution Si B to the action (3.1) to be 



= ^ [ d 3 x $f <P A e AB , (3.8) 

1 JdM 



where h = det(hij) and hij gives the intrinsic 3-metric of the boundary = 1,2,3) , 
while eAB is the alternating spinor [22]. One can then verify that the total action Si is 

extremised at a classical solution, given a specification of $f at the boundary. These are 
our local boundary conditions. Note their similarity to the spin-1 (Maxwell) boundary 
conditions, involving the data = —2in B B , e AB k P>k on E/ and T,p in Eq.(6.18) of 

[8] , which correspond to fixing the (spatial) magnetic-field components P>k ■ Again, for 
spin-2 (gravitational-wave) perturbations, the boundary conditions correspond to fixing a 
projection of the Weyl tensor - see Sec. 10 of [9]. 

At a classical solution, by virtue of the Weyl equations (3.4,5), the volume integral 
vanishes in Eq.(3.1) for the action . Hence, the classical spin-i part of the action resides 
on the boundary and, at lowest (quadratic) order in fermions, is given by 

(^3) class (^2^) class ' (^•"^) 

where Si B is given in Eq.(3.8). 

For integer-spin perturbations, as in [5,6,8,9], we chose vanishing perturbative initial 
data on the initial space-like hypersurface £/ , to simplify the exposition. Correspondingly, 
for massless spin-| perturbations, we take the homogeneous (Dirichlet-type) boundary 
conditions at E/ to be 

$> A (x) = , (3.10) 
(choosing e = — 1 here). Equivalently, 

<j) A = V2 n A A , 4> A ' (3.11) 

on £/ . Hence, there is no contribution to the classical action from the initial surface. 
On the final hypersurface £f, we specify (non-zero) data ^ A (x) , corresponding to the 
presence of a remnant flux of neutrinos outgoing at large radius and late time: 

= f A (r,n) . (3.12) 
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In this case, unless one has found a complete solution of the field equations, ^ is (as yet) 
unknown on Ep, while is unknown on Ej . 

Following the discussion in Sec. 2, the particular choice (3.10), with zero right-hand 
side, would be interpreted as containing no particles on Ej . But the choice (3.12) of final 
data, for f A ^ , specifies that a certain amount of fermions were produced, to arrive 
at T,p . Of course, one can specify the right-hand sides of Eqs.(3.10,12) at will (subject 
to fall-off at infinity) to investigate the effects of allowing incoming fermions on Ej and 
asking how they contribute to the final fermionic output at E^ . Because the classical 
fermionic action in Eq.(4.35) below is 'bosonic' - that is, it is analogous to the particular 
wave function /(a*, a) = fn aa* in Eq.(2.2) - fermions can only be produced or destroyed 
in pairs. For example, from Eq.(4.35) it is impossible to have precisely one neutrino in the 
far future if none were present in the far past. 

In our case, fermionic boundary data are to be specified on a pair of asymptotically- 
flat space-like hypersurfaces E/ and E p , together with the boundary data, which include 
the asymptotic proper-time interval T, as measured at spatial infinity. It turns out [11] 
that one should then specify ^ on one surface and on the other, to have a well-posed 
classical problem. In the 'cosmological' case, where, for example, a three-sphere S 3 bounds 
a four-ball B 4 , different considerations may apply [11,29-31]. 

4. Summary of the Vaidya approximation 

By analogy with the s = case in [1-7] , we now study weak s = \ massless (neutrino) 
fields in a spherically-symmetric 'background' Lorentzian-signature geometry 



along the lines described in [1,2]. We are principally interested in quantum emission of a 
roughly stochastic type, following gravitational collapse to a black hole, since generically 
this leads to the most probable final configurations. At late times, as described further 
in [1-9], the configuration of all fields (not just the 'background' spherically-symmetric 
coupled solution of the Einstein/massless-scalar field equations, but also the averaged 
second-variation energy- momentum tensor T^ v generated by the anisotropic perturbative 
fields) will approximate closely a nearly-spherically-symmetric stream of radially-outgoing 
null or massless radiation, together with the gravitational field which this generates. This 
is treated in detail in [3,4]. Here, we summarise briefly the necessary properties of such 
an approximately Vaidya metric [7,27]. In the companion paper [4] on spins 1 and 2, a 
full description of the aspects of the Vaidya solution relevant to the present spin-1/2 case 
is given. In [7], moreover, a detailed treatment of the Vaidya metric, as needed for the 
general question of black-hole evaporation, can be found. In addition to treatment of the 
Vaidya-like geometry in which the particles or fields propagate, [7] also includes a thorough 
description of the energy-momentum source for the Vaidya-like geometry, produced by a 
large number of null (radiating) particles. 

As in [1-9], we write e~ a ^^ in terms of the 'mass function' m(t, r) as 



ds * = _ e b(t,r) df 2 + e a(t,r) ^2 + r 2 ^2 + ^2 Q # 2 j 



(4.1) 



e 



a(t,r) 



= 1 - 



2m(t, r) 



(4.2) 



r 
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A Vaidya metric [7,27] corresponds to a spherically-symmetric null-fluid energy-momentum 
source, of the form 

Tfiv = /(t,r)fc M fc v ■ (4.3) 

Here, f(t,r) is some spherically-symmetric function, and the vector field /c M is null 
(k^kfj, = 0) and points radially outward to the future (in our case). As is described much 
more fully in the companion paper [4], the corresponding Vaidya solution of the Einstein 
field equations with the energy-momentum source (4.3) can be put in the diagonal form 
(4.1), as 



ds 2 = 



(^) 2 (l-^^) _1 ^ 2 + (l-H^) _1 dr 2 + r 2 dfi 2 . (4.4) 



Here, m = m(t,r) , and m, m' denote dm/dt, dm/dr , respectively. Thus, 



Moreover, as seen in [4,7], the Einstein field equations imply that 

m' (l " ^) = /M , (4-6) 

where /(m) > is an arbitrary (smooth) function. Physically, in our example, the function 
f(m) in Eq.(4.6) depends on the specific physical (particle) model. This metric describes 
the space-time region in which radiation is streaming radially outwards, due to the black- 
hole evaporation. Correspondingly, the 'mass function' m(t, r) varies extremely slowly with 
time, and effectively one has a Schwarzschild metric with a very slowly-varying mass, in 
which f(m) can be regarded as nearly constant. As described in the preceding paper [4], 
the geometry will gradually deviate from the Vaidya form as one moves to the past into 
the strongly dynamical collapse region. Provided that the complexified classical boundary- 
value problem is well-posed, the complexified classical solution (<7^,0) will be regular at 
every point, including the spatial origin r = . 

As further described in [4,7], a coordinate transformation can be found, which leads 
to the usual form 

ds 2 = - ( 1 - \du 2 - 2dudr + r 2 {d9 2 + sin 2 0# 2 ) , (4.7) 



of the Vaidya metric [7,27], where m = m{u) is now seen to be a function of a single 
null variable u. The diagonal form (4.4) will be used below, in studying the evolution of 
linearised perturbation fields on the slowly-varying sequence of Schwarzschild geometries 
which is described closely here by a Vaidya metric. The 'null' form (4.7) is instead useful 
in understanding the large-scale structure of the Vaidya space-time. 



9 



From the Vaidya metric in the form (4.7), one can see that the apparent singularity 
in the metric (4.4) at r = 2m(u) is only a coordinate singularity [27]. Further, the surface 
{r = 2m(u)} is space-like, lying to the past of the region {r > 2m(«)}. In fact, the 
geometry in the region {r < 2m(u)} will gradually deviate from the Vaidya form, as one 
moves to the past by (say) reducing u while holding r fixed, since one reaches the region of 
strong-field gravitational collapse. This region can still be described by the diagonal metric 
(4.1), with scalar field <p = <p(t, r) , but the full field equations enforce a more complicated 
coupled solution. Again, provided that the complexified boundary- value problem of Sec.l 
remains well-posed for < 6 < tv/2 , the full, complexified Einstein/scalar classical solution 
studied here will be regular at the spatial origin r = . 

5. Separation of neutrino equations 

Separation of variables for the neutrino equation in a Kerr background space-time was 
carried out by Teukolsky [32] and by Unruh [33], using the Newman-Penrose formalism 
[34]. As in [5,6] for s = 0, [8] for s = 1 and [9] for s = 2 , we are principally concerned 
here with the wave-like (Lorentzian) evolution of neutrino fields in the Vaidya-like metric 
[7,27], which describes the extremely slow evolution of the black hole, as it gradually loses 
mass through nearly-isotropic radiation. For most radiation frequencies of interest, the 
time rate of change of the background geometry will be so slow as to be negligible. Thus, 
the evolution of the wave will be adiabatic. 

Following [32,33], we decompose the neutrino field (4> A , <p A ) with respect to the nor- 
malised spinor dyad (t A , o A ) which corresponds to the Kinnersley null tetrad [35] in the 
Schwarzschild metric. This pair (o A , i A ) gives a basis for the 2-complex-dimensional vec- 
tor space of spinor s u A at each point, and is normalised according to oa t A = 1 = 
— la o a . Knowledge of the Kinnersley null tetrad [34,35] , n M , , of vectors 
at a point is equivalent to knowledge of the corresponding normalised spinor dyad (o A , l a ). 




(* + ih d * + \ cot ") xW + \ ( a < ~ e ~° 9 - ~ i - L 7^ l ) x ' 1 ' = • <6 ' 1) 

H< + a - + ;) * <0) + ^ ( de + ^> a * + \ cote ) xW = - <6 ' 2) 



As is standard, one can separate variables in Eqs.(5.1,2). We write 



X 




(5.3) 



X 




(5.4) 
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From Eqs. (5.3,4), S ± i tm (6,<f>) and R±i £m (9,(j)) satisfy: 



1 .2 ( a „-« a (r-M) 



-r ( d t - e a d r - 2 — J R+i em = - Xi R-^ m , (5.5) 



If ml \ 

7f I 96 ~ ^9 + 2 COte ) S ~¥m = A l S +$lm > ( 5 - 6 ) 

(e a d t + d r )R_ hem = -X 2 R +htn , (5.7) 

If ml \ 

^{ d9+ sh[9 + 2 C ° te ) S+1 ^ m = - X * S -¥™ > ( 5 - 8 ) 

where Ai and A2 are separation constants. Regularity of the angular functions S±i im at 
9 = and 9 = n implies that Ai = A2 , = Xi , say. That is, if we set 6' = (tt — 9) in 
Eq.(5.6,8), we find that S±i im (9) oc iS l = pi^ m (7r — 9) . Up to normalisation, the functions 

S±i lm e im< ^ are just spin-weighted spherical harmonics s >£ m (fi) for s = ±| [36,37]. 
These functions have the properties 

.Ytm = (-!) m+S -Yt,-m , (5-9) 

s Y em (n-0, 7T + 0) = s ^ m (0,0) , (5.10) 

y dfi 8 l£ m (f2) fl l^ m ,(f2) = 8 U > 5 mm > . (5.11) 
In general, the S s £ m , which are regular on [0, it] , satisfy: 



1 d f . dS 8 i m 
siny 



sin# d0V <20 



/ — m 2 2ms cot #99^ \ „ 

+ — — n — Q S 2 COt 2 9 + S + V) S sim = , 

V sin 9 sin 9 J 

where v = (£ — s)(£ + s + 1) . For our case \s\ = \ , one has v = 2 (A^) 2 = (£ + 
Finally, we make a further change of variables to 



(5.12) 



1)2 



and then to 



P 



P+±Ut,r) = ^=e~ a / 2 R Hem (t,r) , (5.13) 
P-Um& r ) = R -Um(t,r) , (5.14) 



Um(t,r) = 1== Ui tm +(t,r) - Zu m -(t,rj) , (5.15) 

P+Um(t,r) = I U^ m +(t,r) + £i em _(t,rj) . (5.16) 
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We obtain the more symmetrical set of equations 



e - Q r + -^jJL e~ a ^ £i, m _ = d t £ 



e~ a d r - '-^ e""/ 2 ) ti tm+ = d t i ¥m _ . 

These yield the second-order decoupled (t, r) form of the neutrino wave equation 

e~ a d r (e~ a d r ^i em± ) - (d t f fi £m± - Vi e ± £i em ± = , 
where the real potential(s) Vi e ± are given by 



(5.17) 
(5.18) 

(5.19) 



Vi 



1/2 



(5.20) 



Thus, we have expressed the perturbative neutrino field in terms of 'odd' and 'even', 
or 'magnetic' and 'electric' components £i^ m ± , just as we did with the integer-spin per- 
turbations. One can now use the Kinnersley tetrad 



r = (1/(1-^), 1,0,0) , 
- (i,-(i-^),o,o)/2 , 

n? = ^0,0,1, (i/sm8yj/rV2 , 



(5.21) 



in the Schwarzschild metrics together with the boundary condition (3.10), to write the 
classical action (3.8,9) as 



(Si s)class 



dQ 

'o 

00 £ 



Roo 



dr r 2 e a ' 2 o B X (0) + i B 



= 4j E E 1^ dr e* i Um _ t iim . 



t— 2 

Here we have also used Eqs.(5.9,ll,15,16). 

The boundary condition (3.10) translates into 



(5.22) 



From Eq.(5.21), one finds 



= . 



t=o 



= . 



t=0 
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(5.23) 



(5.24) 



In addition, one has 

Hoa = + V2e^x^ = £ H Y ^) ,(5.25) 



* Ivn \/£ + o 



2 

Accordingly, we set 

^i^ m _(t,r) = j & ci Mm _ ^i u _(r) ^ , (5.27) 

where the {ci Wm _} are certain coefficients taking (odd) values in a Grassmann algebra, 
and the {£ik£-( r )} are radial functions. As above, we have assumed that, in a neighbour- 
hood of Hp , an adiabatic approximation is valid, such that the frequencies are large when 
compared with the inverse time-scale for the very slow change in the background metric. 

From Eqs.(5. 17-19), regularity near the centre of symmetry {r = 0} demands that 

Z$kie(r) ~ Vkee r j e _i e (kr) , (5.28) 

as r — > , where ji_i e (kr) denotes a spherical Bessel function [38], and fx k £e is a normal- 
isation constant. Hence, the {£i^ e (?")} are rea ^ functions. 
Noting the (£, r) neutrino field equation (5.19), we set 

^u-ir) = ({*lu-)e ikr * + (*!«-) e- ifcr *)/i M _(r) , (5-29) 

where the zi kt _ are complex coefficients, and we take the boundary condition at infinity 

to be /i M _(r) -> 1 as r -> oo , so that Ci^-W ~ ((^/w-) ^ + e" lfcr *) 
as r — > oo . Here, r* is a generalised Regge- Wheeler coordinate [39,40], such that 

^ = e^/ 2 , (5.30) 

and one also has fi^_ ki _(r) = fi k£ _(r). The coefficients <zi fc ^_ satisfy z\ u _ = 
zi t -k£- • The normalisation of the radial functions is given by 

J q R °° dr e a fi M _(r) Ci^-W = 2tt zi m _ (*(*, k') + 8(k, -A/)) . (5.31) 

Hence, the classical spin-i action (5.22) can be written in the form 



(Si B ) c i ass = -inV2j2 £ 



x 



X 



(5.32) 
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This is the desired form of the spin--| classical action, providing the nearest corre- 
spondence with the bosonic expressions Eq.(3.8) of [6], Eq.(5.19) of [10] and Eq.(9.7) of 
[9] for spins s = , 1 and 2 , respectively. When one compares this with the massless- 
scalar (second- variation) action S c \ a ss in Eq.(3.8) of [6], the s = expression involves an 
integral f_ dk k ( ) , whereas for s = ^ we have J °° dw ( ) , with no power of u . The 
s = integrand is schematically of the form zzaa , whereas for s = \ one has z zee, 
where the tildes refer to the primed spinor field 4> A ' , and untilded quantities refer to the 
unprimed field 4> A ; this structure is inevitable in massless fermionic theories. Finally, the 
s = integrand includes a factor cot(/cT), which was crucial in determining the form of 
the s — quantum amplitude in [6] , via rotation of T into the complex. In the s = ^ case, 
this factor of cot(fcT) is absent, and the calculation of the s = | quantum amplitude is 
correspondingly simpler. 

6. Conclusion 

We have seen in this paper how the calculations in [5-9], treating quantum ampli- 
tudes for the emission of bosonic particles (s = 0, 1, 2), in a nearly-spherically-symmetric 
Einstein/massless-scalar collapse to a black hole, can be extended to the fermionic case of 
a massless spin-i neutrino field. Following the treatment in [1-3,5-7] of the s = massless- 
scalar case, the classical action functional for s = \ data on the final late-time space-like 
hypersurface Ep has been calculated. One again takes the initial Einstein/scalar data on 
E/ to be exactly spherically symmetric, describing (say) diffuse, slowly-moving initial data 
for gravitational collapse to a black hole. We assume that the (real) Lorentzian proper-time 
separation T between E/ and Ep , as measured at spatial infinity, is sufficiently large that 
Ei? intersects all the radiation emitted. Moreover, the final data consist of a section of 
a spherically-symmetric approximately Vaidya geometry, together with weak spin--| data. 
We rotate T — >• \T\ exp(—i6) into the lower-half complex-T plane (0 < 9 < 7r/2), after Feyn- 
man. We expect that the classical boundary-value problem will then be well-posed. The 
classical action has here been calculated as a functional of the (suitably chosen) boundary 
data, leading to the quantum amplitude. 

Perhaps the main differences between the bosonic calculations and the present fermionic| 
calculation are that, classically, the spin-i field is described by odd (anti-commuting) 
Grassmann quantities, and, of course, that the field equations are of first order, unlike the 
commuting and second-order nature of the bosonic systems. The first-order property for 
spin-i is reflected in our choice of local boundary conditions, 'inherited' from local su- 
persymmetry [24-26,41,42]. This 'ancestral' appearance of supergravity is reflected in the 
form of the classical action functional of the corresponding 'local' boundary data, as well 
as in the resulting quantum amplitude. An analogous (but more complicated) treatment 
of amplitudes for the spin-| field, crucial for local supersymmetry (supergravity) will soon 
be completed [43]. 

As remarked in [2], the present approach leads to a calculation of quantum amplitudes, 
appropriate for studying the black-hole radiation which follows gravitational collapse. This 
formulation must be different from the familiar one (which is normally carried out by 
considering Bogoliubov transformations), since it yields quantum amplitudes relating to 
the final state, and not just the usual probabilities for outcomes at a late time and large 
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radius. It is not that many of the attributes of the radiation, whether described in the 
present way or in the usual way, can be very different in the two descriptions. It is rather 
that the present description allows one to ask more questions of the radiation (and to 
receive a response) . 
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